A Solanum lycopersicoides reference genome facilitates insights into tomato specialized metabolism and immunity.

Wild relatives of tomato are a valuable source of natural variation in tomato breeding, as many can be hybridized to the cultivated species (Solanum lycopersicum). Several, including Solanum lycopersicoides, have been crossed to S. lycopersicum for the development of ordered introgression lines (ILs), facilitating breeding for desirable traits. Despite the utility of these wild relatives and their associated ILs, few finished genome sequences have been produced to aid genetic and genomic studies. Here we report a chromosome-scale genome assembly for S. lycopersicoides LA2951, which contains 37 938 predicted protein-coding genes. With the aid of this genome assembly, we have precisely delimited the boundaries of the S. lycopersicoides introgressions in a set of S. lycopersicum cv. VF36 × LA2951 ILs. We demonstrate the usefulness of the LA2951 genome by identifying several quantitative trait loci for phenolics and carotenoids, including underlying candidate genes, and by investigating the genome organization and immunity-associated function of the clustered Pto gene family. In addition, syntenic analysis of R2R3MYB genes sheds light on the identity of the Aubergine locus underlying anthocyanin production. The genome sequence and IL map provide valuable resources for studying fruit nutrient/quality traits, pathogen resistance, and environmental stress tolerance. We present a new genome resource for the wild species S. lycopersicoides, which we use to shed light on the Aubergine locus responsible for anthocyanin production. We also provide IL boundary mappings, which facilitated identifying novel carotenoid quantitative trait loci of which one was likely driven by an uncharacterized lycopene ß-cyclase whose function we demonstrate.

Whole genome sequencing and assembly data of Moricandia moricandioides and M. arvensis.

Moricandia is a genus belonging to the family Brassicaceae. C3 and C3-C4 photosynthesis Moricandia species exist in a close phylogenetic proximity, as well as to Brassica crops. Here, we performed PacBio genome sequencing on M. moricandioides and M. arvensis. The genomes were assembled using Flye assembler, then polished with Illumina reads and reduced duplication with Purge Haplotigs. The total length of genome assemblies of M. moricandioides and M. arvensis was 498 Mbp and 759 Mbp, respectively. These data will be useful for studies of the genetic control of C3-C4 characteristics, therefore gaining new insights into the early evolutionary steps of C4 photosynthesis. Further, it can be integrated into Brassica crop breeding. The data can be accessed at ENA under the project number PRJEB39764.

Helixer: cross-species gene annotation of large eukaryotic genomes using deep learning.

MOTIVATION: Current state-of-the-art tools for the de novo annotation of genes in eukaryotic genomes have to be specifically fitted for each species and still often produce annotations that can be improved much further. The fundamental algorithmic architecture for these tools has remained largely unchanged for about two decades, limiting learning capabilities. Here, we set out to improve the cross-species annotation of genes from DNA sequence alone with the help of deep learning. The goal is to eliminate the dependency on a closely related gene model while also improving the predictive quality in general with a fundamentally new architecture. RESULTS: We present Helixer, a framework for the development and usage of a cross-species deep learning model that improves significantly on performance and generalizability when compared to more traditional methods. We evaluate our approach by building a single vertebrate model for the base-wise annotation of 186 animal genomes and a separate land plant model for 51 plant genomes. Our predictions are shown to be much less sensitive to the length of the genome than those of a current state-of-the-art tool. We also present two novel post-processing techniques that each worked to further strengthen our annotations and show in-depth results of an RNA-Seq based comparison of our predictions. Our method does not yet produce comprehensive gene models but rather outputs base pair wise probabilities. AVAILABILITY AND IMPLEMENTATION: The source code of this work is available at https://github.com/weberlab-hhu/Helixer under the GNU General Public License v3.0. The trained models are available at https://doi.org/10.5281/zenodo.3974409. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Phosphate starvation causes different stress responses in the lipid metabolism of tomato leaves and roots.

Plants have evolved various acclimation responses to cope with phosphate depletion, including several changes in lipid metabolism. Thereby membrane phospholipids are dephosphorylated and can be used as an internal phosphate source, while galactolipids are incorporated into the membrane to maintain membrane functionality. Still little is known about the lipidomic and transcriptomic response of plants other than Arabidopsis thaliana upon phosphate starvation. Therefore, we employed lipidomics and transcriptomics to characterize the phosphate starvation response of lipid metabolism in tomato leaves and roots. Overall, phospholipid levels decreased and galactolipids increased during the acclimation response. In addition, an early increase of triacylglycerol was observed. Interestingly, there were major differences in the acclimation response of tomato leaves and roots: leaves mainly accumulated polyunsaturated triacylglycerol, while roots showed a massive increase in galactolipid content. In line with these results, we observed transcriptional induction of phospholipid degradation and galactolipid synthesis pathways in both analyzed tissues. In contrast, other aspects of the transcriptional response, in particular, the induction of phospholipid degradation, ER-localized fatty acid desaturation and triacylglycerol assembly differed between tomato leaves and roots. These results suggest a different modulation of degraded phospholipids toward triacylglycerols and galactolipids in phosphate-starved tomato leaves and roots. Possibly the availability and composition of acyl-CoA pools and ER-derived precursors trigger the synthesis of triacylglycerols or galactolipids. As the mechanism of triacylglycerol accumulation is poorly characterized outside of seed oil formation, these findings enhance our understanding of the phosphate starvation response and of how storage lipids accumulate under stress in vegetative tissue.

Transcriptomic and proteomic data in developing tomato fruit.

Transcriptomic and proteomic analyses were performed on three replicates of tomato fruit pericarp samples collected at nine developmental stages, each replicate resulting from the pooling of at least 15 fruits. For transcriptome analysis, Illumina-sequenced libraries were mapped on the tomato genome with the aim to obtain absolute quantification of mRNA abundance. To achieve this, spikes were added at the beginning of the RNA extraction procedure. From 34,725 possible transcripts identified in the tomato, 22,877 were quantified in at least one of the nine developmental stages. For the proteome analysis, label-free liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) was used. Peptide ions, and subsequently the proteins from which they were derived, were quantified by integrating the signal intensities obtained from extracted ion currents (XIC) with the MassChroQ software. Absolute concentrations of individual proteins were estimated for 2375 proteins by using a mixed effects model from log10-transformed intensities and normalized to the total protein content. Transcriptomics data are available via GEO repository with accession number GSE128739. The raw MS output files and identification data were deposited on-line using the PROTICdb database (http://moulon.inra.fr/protic/tomato_fruit_development) and MS proteomics data have also been deposited to the ProteomeXchange with the dataset identifier PXD012877. The main added value of these quantitative datasets is their use in a mathematical model to estimate protein turnover in developing tomato fruit.

Modeling Protein Destiny in Developing Fruit.

Protein synthesis and degradation are essential processes that regulate cell status. Because labeling in bulky organs, such as fruits, is difficult, we developed a modeling approach to study protein turnover at the global scale in developing tomato (Solanum lycopersicum) fruit. Quantitative data were collected for transcripts and proteins during fruit development. Clustering analysis showed smaller changes in protein abundance compared to mRNA abundance. Furthermore, protein and transcript abundance were poorly correlated, and the coefficient of correlation decreased during fruit development and ripening, with transcript levels decreasing more than protein levels. A mathematical model with one ordinary differential equation was used to estimate translation (kt ) and degradation (kd ) rate constants for almost 2,400 detected transcript-protein pairs and was satisfactorily fitted for >1,000 pairs. The model predicted median values of ∼2 min for the translation of a protein, and a protein lifetime of ∼11 d. The constants were validated and inspected for biological relevance. Proteins involved in protein synthesis had higher kt and kd values, indicating that the protein machinery is particularly flexible. Our model also predicts that protein concentration is more strongly affected by the rate of translation than that of degradation.

Footprints of parasitism in the genome of the parasitic flowering plant Cuscuta campestris.

A parasitic lifestyle, where plants procure some or all of their nutrients from other living plants, has evolved independently in many dicotyledonous plant families and is a major threat for agriculture globally. Nevertheless, no genome sequence of a parasitic plant has been reported to date. Here we describe the genome sequence of the parasitic field dodder, Cuscuta campestris. The genome contains signatures of a fairly recent whole-genome duplication and lacks genes for pathways superfluous to a parasitic lifestyle. Specifically, genes needed for high photosynthetic activity are lost, explaining the low photosynthesis rates displayed by the parasite. Moreover, several genes involved in nutrient uptake processes from the soil are lost. On the other hand, evidence for horizontal gene transfer by way of genomic DNA integration from the parasite's hosts is found. We conclude that the parasitic lifestyle has left characteristic footprints in the C. campestris genome.

De Novo Assembly of a New Solanum pennellii Accession Using Nanopore Sequencing.

Updates in nanopore technology have made it possible to obtain gigabases of sequence data. Prior to this, nanopore sequencing technology was mainly used to analyze microbial samples. Here, we describe the generation of a comprehensive nanopore sequencing data set with a median read length of 11,979 bp for a self-compatible accession of the wild tomato species Solanum pennellii We describe the assembly of its genome to a contig N50 of 2.5 MB. The assembly pipeline comprised initial read correction with Canu and assembly with SMARTdenovo. The resulting raw nanopore-based de novo genome is structurally highly similar to that of the reference S. pennellii LA716 accession but has a high error rate and was rich in homopolymer deletions. After polishing the assembly with Illumina reads, we obtained an error rate of <0.02% when assessed versus the same Illumina data. We obtained a gene completeness of 96.53%, slightly surpassing that of the reference S. pennellii Taken together, our data indicate that such long read sequencing data can be used to affordably sequence and assemble gigabase-sized plant genomes.

Early metabolic and transcriptional variations in fruit of natural white-fruited Fragaria vesca genotypes.

Strawberry fruits (Fragaria vesca) are valued for their sweet fruity flavor, juicy texture, and characteristic red color caused by anthocyanin pigments. To gain a deeper insight into the regulation of anthocyanin biosynthesis, we performed comparative metabolite profiling and transcriptome analyses of one red-fruited and two natural white-fruited strawberry varieties in two tissues and three ripening stages. Developing fruit of the three genotypes showed a distinctive pattern of polyphenol accumulation already in green receptacle and achenes. Global analysis of the transcriptomes revealed that the ripening process in the white-fruited varieties is already affected at an early developmental stage. Key polyphenol genes showed considerably lower transcript levels in the receptacle and achenes of both white genotypes, compared to the red genotype. The expression of the anthocyanidin glucosyltransferase gene and a glutathione S-transferase, putatively involved in the vacuolar transport of the anthocyanins, seemed to be critical for anthocyanin formation. A bHLH transcription factor is among the differentially expressed genes as well. Furthermore, genes associated with flavor formation and fruit softening appear to be coordinately regulated and seem to interact with the polyphenol biosynthesis pathway. This study provides new information about polyphenol biosynthesis regulators in strawberry, and reveals genes unknown to affect anthocyanin formation.

Freeze-quenched maize mesophyll and bundle sheath separation uncovers bias in previous tissue-specific RNA-Seq data.

The high efficiency of C4 photosynthesis relies on spatial division of labor, classically with initial carbon fixation in the mesophyll and carbon reduction in the bundle sheath. By employing grinding and serial filtration over liquid nitrogen, we enriched C4 tissues along a developing leaf gradient. This method treats both C4 tissues in an integrity-preserving and consistent manner, while allowing complementary measurements of metabolite abundance and enzyme activity, thus providing a comprehensive data set. Meta-analysis of this and the previous studies highlights the strengths and weaknesses of different C4 tissue separation techniques. While the method reported here achieves the least enrichment, it is the only one that shows neither strong 3' (degradation) bias, nor different severity of 3' bias between samples. The meta-analysis highlighted previously unappreciated observations, such as an accumulation of evidence that aspartate aminotransferase is more mesophyll specific than expected from the current NADP-ME C4 cycle model, and a shift in enrichment of protein synthesis genes from bundle sheath to mesophyll during development. The full comparative dataset is available for download, and a web visualization tool (available at http://www.plant-biochemistry.hhu.de/resources.html) facilitates comparison of the the Z. mays bundle sheath and mesophyll studies, their consistencies and their conflicts.

Understanding metabolite transport and metabolism in C4 plants through RNA-seq.

RNA-seq, the measurement of steady-state RNA levels by next generation sequencing, has enabled quantitative transcriptome analyses of complex traits in many species without requiring the parallel sequencing of their genomes. The complex trait of C4 photosynthesis, which increases photosynthetic efficiency via a biochemical pump that concentrates CO2 around RubisCO, has evolved convergently multiple times. Due to these interesting properties, C4 photosynthesis has been analyzed in a series of comparative RNA-seq projects. These projects compared both species with and without the C4 trait and different tissues or organs within a C4 plant. The RNA-seq studies were evaluated by comparing to earlier single gene studies. The studies confirmed the marked changes expected for C4 signature genes, but also revealed numerous new players in C4 metabolism showing that the C4 cycle is more complex than previously thought, and suggesting modes of integration into the underlying C3 metabolism.

Retained duplicate genes in green alga Chlamydomonas reinhardtii tend to be stress responsive and experience frequent response gains.

BACKGROUND: Green algae belong to a group of photosynthetic organisms that occupy diverse habitats, are closely related to land plants, and have been studied as sources of food and biofuel. Although multiple green algal genomes are available, a global comparative study of algal gene families has not been carried out. To investigate how gene families and gene expression have evolved, particularly in the context of stress response that have been shown to correlate with gene family expansion in multiple eukaryotes, we characterized the expansion patterns of gene families in nine green algal species, and examined evolution of stress response among gene duplicates in Chlamydomonas reinhardtii. RESULTS: Substantial variation in domain family sizes exists among green algal species. Lineage-specific expansion of families occurred throughout the green algal lineage but inferred gene losses occurred more often than gene gains, suggesting a continuous reduction of algal gene repertoire. Retained duplicates tend to be involved in stress response, similar to land plant species. However, stress responsive genes tend to be pseudogenized as well. When comparing ancestral and extant gene stress response state, we found that response gains occur in 13% of duplicate gene branches, much higher than 6% in Arabidopsis thaliana. CONCLUSION: The frequent gains of stress response among green algal duplicates potentially reflect a high rate of innovation, resulting in a species-specific gene repertoire that contributed to adaptive response to stress. This could be further explored towards deciphering the mechanism of stress response, and identifying suitable green algal species for oil production.

Comparative transcriptome atlases reveal altered gene expression modules between two Cleomaceae C3 and C4 plant species.

C(4) photosynthesis outperforms the ancestral C(3) state in a wide range of natural and agro-ecosystems by affording higher water-use and nitrogen-use efficiencies. It therefore represents a prime target for engineering novel, high-yielding crops by introducing the trait into C(3) backgrounds. However, the genetic architecture of C(4) photosynthesis remains largely unknown. To define the divergence in gene expression modules between C(3) and C(4) photosynthesis during leaf ontogeny, we generated comprehensive transcriptome atlases of two Cleomaceae species, Gynandropsis gynandra (C(4)) and Tarenaya hassleriana (C(3)), by RNA sequencing. Overall, the gene expression profiles appear remarkably similar between the C(3) and C(4) species. We found that known C(4) genes were recruited to photosynthesis from different expression domains in C(3), including typical housekeeping gene expression patterns in various tissues as well as individual heterotrophic tissues. Furthermore, we identified a structure-related module recruited from the C(3) root. Comparison of gene expression patterns with anatomy during leaf ontogeny provided insight into genetic features of Kranz anatomy. Altered expression of developmental factors and cell cycle genes is associated with a higher degree of endoreduplication in enlarged C(4) bundle sheath cells. A delay in mesophyll differentiation apparent both in the leaf anatomy and the transcriptome allows for extended vein formation in the C(4) leaf.

Predicting C4 photosynthesis evolution: modular, individually adaptive steps on a Mount Fuji fitness landscape.

An ultimate goal of evolutionary biology is the prediction and experimental verification of adaptive trajectories on macroevolutionary timescales. This aim has rarely been achieved for complex biological systems, as models usually lack clear correlates of organismal fitness. Here, we simulate the fitness landscape connecting two carbon fixation systems: C3 photosynthesis, used by most plant species, and the C4 system, which is more efficient at ambient CO2 levels and elevated temperatures and which repeatedly evolved from C3. Despite extensive sign epistasis, C4 photosynthesis is evolutionarily accessible through individually adaptive steps from any intermediate state. Simulations show that biochemical subtraits evolve in modules; the order and constitution of modules confirm and extend previous hypotheses based on species comparisons. Plant-species-designated C3-C4 intermediates lie on predicted evolutionary trajectories, indicating that they indeed represent transitory states. Contrary to expectations, we find no slowdown of adaptation and no diminishing fitness gains along evolutionary trajectories.

C4 photosynthesis: from evolutionary analyses to strategies for synthetic reconstruction of the trait.

C4 photosynthesis represents the most productive modes of photosynthesis in land plants and some of the most productive crops on the planet, such as maize and sugarcane, and many ecologically important native plants use this type of photosynthesis. Despite its ecological and economic importance, the genetic basis of C4 photosynthesis remains largely unknown. Even many fundamental aspects of C4 biochemistry, such as the molecular identity of solute transporters, and many aspects of C4 plant leaf development, such as the Kranz anatomy, are currently not understood. Here, we review recent progress in gaining a mechanistic understanding of the complex C4 trait through comparative evolutionary analyses of C3 and C4 species.

Systems analysis of a maize leaf developmental gradient redefines the current C4 model and provides candidates for regulation.

We systematically analyzed a developmental gradient of the third maize (Zea mays) leaf from the point of emergence into the light to the tip in 10 continuous leaf slices to study organ development and physiological and biochemical functions. Transcriptome analysis, oxygen sensitivity of photosynthesis, and photosynthetic rate measurements showed that the maize leaf undergoes a sink-to-source transition without an intermediate phase of C(3) photosynthesis or operation of a photorespiratory carbon pump. Metabolome and transcriptome analysis, chlorophyll and protein measurements, as well as dry weight determination, showed continuous gradients for all analyzed items. The absence of binary on-off switches and regulons pointed to a morphogradient along the leaf as the determining factor of developmental stage. Analysis of transcription factors for differential expression along the leaf gradient defined a list of putative regulators orchestrating the sink-to-source transition and establishment of C(4) photosynthesis. Finally, transcriptome and metabolome analysis, as well as enzyme activity measurements, and absolute quantification of selected metabolites revised the current model of maize C(4) photosynthesis. All data sets are included within the publication to serve as a resource for maize leaf systems biology.